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CP-PACS and JLQCD Collaborations are carrying out a joint project of the 2-1-1 flavor full QCD with the RG- 
improved gauge action and the non-perturbatively C'(a)-improved Wilson quark action. This simulation removes 
quenching effects of all three light quarks, which is the last major uncertainty in lattice QCD. In this report we 
present our results for the light meson spectrum and quark masses on a 20^ x 40 lattice at the lattice spacing 
a ~ 0.10 fm. 
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1. Introduction 

Three-flavor full QCD simulation is a current 
crucial task of lattice QCD, as it would enable 
us, for the first time, to study hadron physics 
based on principles of QCD with no approxima- 
tions. To achieve this, the CP-PACS and JLQCD 
Collaborations are carrying out simulations with 
dynamical degenerate u, d quarks and a s quark 
with different quark mass. 

In our previous work presented at Lattice 2003, 
we made a exploratory study at a lattice spacing 
a ~ 0.10 fm on a small lattice 16'^ x 32 [1]. In this 
article we report the new results for the light me- 
son spectrum and quark masses on a larger lattice 
of 20^ X 40 at a ~ 0.10 fm. We note that the anal- 
ysis procedures have been made more systematic 
over these adopted at the time of Lattice 2004, 
and results are, albeit qualitatively the same, nu- 
merically different. 
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2. Strategy of the project 

We employ the RG-improved gauge action and 
the non-perturbatively 0{a) improved Wilson 
quark action determined for this gauge action [2] , 
and plan to perform simulations at three lattice 
spacings, a ~ 0.07 fm, 0.10 fm, and 0.122 fm, 
which are at even intervals in a^, with a fixed 
physical volume ^ (2.0fm)'^. We take five u, 
d quark masses for chiral extrapolations in the 
range of mps / my — 0.78 — 0.62, and two s quark 
masses around nips/mv — 0.7 for interpolations. 
These parameters are chosen from our experience 
of simulations at a ~ 0.10 fm on a 16^ X 32 lattice. 

An estimate of computational resources avail- 
able to this project, including the Earth Simu- 
lator, and performance of the polynomial hybrid 
Monte Carlo program [3, 4] we use indicates that 
we can accumulate at least 5000 trajectories at 
a ~ 0.10 fm and 0.122 fm and 2000 trajectories 
at a 2± 0.07 fm, within two to three years, of which 
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Table 1 

At parameters, statistics and measured meson masses (in lattice unit). 





Kg 


traj. 


mps,LL 


mps,LS 




mv,LL 


'mv,LS 




0.13580 


0.13580 


5000 


0.5390(09) 






0.7024(17) 






0.13610 


0.13580 


6000 


0.4943(07) 


0.5098(07) 


0.5249(07) 


0.6656(19) 


0.6764(18) 


0.6871(18) 


0.13640 


0.13580 


7500 


0.4460(06) 


0.4788(06) 


0.5100(05) 


0.6189(19) 


0.6414(17) 


0.6635(15) 


0.13680 


0.13580 


8000 


0.3786(08) 


0.4382(07) 


0.4921(06) 


0.5642(27) 


0.6029(22) 


0.6410(18) 


0.13700 


0.13580 


8000 


0.3384(09) 


0.4140(08) 


0.4799(07) 


0.5300(19) 


0.5776(16) 


0.6243(14) 


0.13580 


0.13640 


5200 


0.5254(09) 


0.4947(10) 


0.4626(11) 


0.6853(14) 


0.6636(15) 


0.6416(17) 


0.13610 


0.13640 


8000 


0.4804(07) 


0.4643(07) 


0.4477(07) 


0.6451(14) 


0.6340(15) 


0.6229(16) 


0.13640 


0.13640 


9000 


0.4324(05) 






0.6045(17) 






0.13680 


0.13640 


9000 


0.3614(06) 


0.3867(06) 


0.4107(06) 


0.5441(18) 


0.5598(16) 


0.5754(15) 


0.13700 


0.13640 


8000 


0.3221(07) 


0.3625(07) 


0.3994(06) 


0.5164(24) 


0.5396(21) 


0.5625(19) 



we are in the second year. 

In this project, our first interest is the meson 
spectrum and light quark masses, because our 
physical volume (2.0fm)'^ is not sufficiently large 
for baryons. 



m,^. (K-input) % (K-input) m,^ (((i-input) m^. (4>-input) 



3. Results at 



0.10 fm 



3.1. Simulations and analyses 

Among three lattice spacings, calculations at 
a ~ 0.10 fm have been completed by simulations 
carried out at /3 = 1.9 with csw — 1.715 on a 
20^^ X 40 lattice. In Tab. ^we show the hopping 
parameters K„d and for u, d (Light) quarks 
and s (Strange) quark. Meson masses are cal- 
culated for Light-Light (LL), Light-Strange (LS) 
and Strange-Strange (SS) combinations of valence 
quark masses, and are given in the table. They 
are determined by a single hyperbolic-cosine cor- 
related fit with tmin — 10. Masses are stable 
against a variation of tmin- 

Chiral fits arc made to LL, LS and SS masses 
simultaneously ignoring correlations among these 
masses. We adopt polynomial functions in quark 
masses up to quadratic terms, x^/d.o.f. (1.1 for 
PS and 0.9 for V) is acceptable. We also test the 
fits using masses normalized by Sommer scale and 
find that the meson spectrum and quark masses 
are stable against the change. 

3.2. Light meson spectrum 

The physical point is determined from ex- 
perimental values of mT^,mp,mK (K-input) or 
mTy,mp,m^ (0-input). The inverse of the lattice 
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Figure 1. Deviation of the meson spectrum from 
experiment at a ~ 0.1 fm. 



spacing = 1.97(4) GeV is independent of in- 
puts for the s quark mass. 

Light meson masses turn out to be 



rriK- = 884.4(3.0) MeV 
= 998.5(5.9) MeV 
rriK = 519.4(6.2) MeV 
ruK' = 895.1(0.3) MeV 



(K-input) , 
(K-input), 
(0- input), 
(0-input). (1) 



In Fig^ we show the deviation of the meson 
spectrum from experiment comparing with our 
results in Nf = 2 and QCD [5] at the same lat- 
tice spacing a 0.1 fm. We find that the spec- 
trum in Nf = 2 + 1 QCD is closer to experiment 
than that in Nf — 2 and QCD at this lattice 
spacing. 
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3.3. Quark masses 

Light quark masses are determined by the 
axial-vector Ward identity. The matching to the 
MS scheme is performed by using the mean-field 
improved one-loop calculation [6] at the scale 
fi = a^^ . The renormalized quark masses are 
evolved to /i = 2 GeV using the 4-loop beta func- 
tion. 

Results for quark masses are 



mud = 3.05(6) 
TO, = 80.4(1.9) 

3.04(6) 
9.3(2.9) 



mud 

TO ,, : 



MeV 
MeV 

MeV 
MeV 



(K-input), 



-input) 



(2) 



Fig. 12 presents quark masses in Nf = 2 + 1 QCD 
comparing with Nf — 2 and QCD at a sim- 
ilar lattice spacing a ^ 0.1 fm. We find that 
quark masses in Nf = 2 + 1 QCD are smaller 
than those in N f = 2 and QCD. Since the non- 
perturbatively improved clover coefficient is used, 
we speculate that the values of quark masses at 
this lattice spacing do not change sizably in the 
continuum limit. 

4. Conclusions and future plans 

As far as we compare the spectrum and quark 
masses at a ~ 0.1 fm, we observe the difference 
between Nf = 2 + 1 and 2. At this stage we 
cannot conclude whether this effect is due to dy- 
namical strange quark or non-perturbatively de- 
termined csvv- This point should become clear in 
the continuum limit. 

Coupling values and hopping parameters for 
the coarser (a ~ 0.122 fm) and finer (a 2± 0.07 fm) 
lattices have already been fixed and production 
runs are in progress. We hope to report results 
in the continuum limit in near future. 

This work is supported by the Epoch Making 
Simulation Projects of Earth Simulator Center, 
the Supercomputer Project No.98 (FY2003) of 
High Energy Accelerator Research Organization 
(KEK), the Large Scale Simulation Projects of 
Academic Computing and Communications Cen- 
ter, University of Tsukuba, Super Sinet Projects 
of National Institute of Infomatics, and also by 
the Grant-in-Aid of the Ministry of Education 
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Figure 2. u and d quark masses (top) and s quark 
mass (bottom) at a 0.1 fm. 
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